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Abstract 

In this work we investigate the influence of substrate temperature on the surface morphology for 
substrate coverage below one monolayer. The model of film growth is based on random deposition 
enriched by limited surface diffusion. Also anisotropy in the growth is involved. We found from 
computer simulations for the simple cubic lattice and solid-on-solid model, that the surface morphology 
changes with increasing temperature from isotropically distributed isolated small islands, through 
anisotropic 1-D stripes to larger 2-D anisotropic islands and again randomly distributed single atoms. 
The transition is also marked in height-height correlation function dependence on temperature, as 
directly seen by snapshots from simulations. The results are in good qualitative agreement with 
already published results of kinetic Monte Carlo simulations, as well as with some experimental data. 

Keywords: Anisotropy; Computer simulation; Epitaxy; Growth; Surface structure, morphology, rough- 
ness and topography. 

1 Introduction 

Surface growth is a problem which may be found in mathematics, physics, biology or even social phenom- 
ena, if we define the surface as a boundary between different phases in d-dimensional space |0, ||, ^, Q . 
Theoretical studies of the properties of the surfaces of the grown crystal/thin films based on Monte Carlo 
(MC) simulations, initially using the solid-on-solid (SOS) models, then extended to more realistic kinetic 
MC methods, became particularly popular in early 80's (see [|| for review). These studies resulted in 
at least a semi-quantitative understanding of growth modes and step dynamics. Then rapid progress in 
observational methods (including diffractional methods — RHEED, LEED and direct imaging — STM, 
AFM) allowed for better understanding mechanism and physics of crystal/films growth. However, MC 
simulations may be still applied as a complementary method which is sometimes more efficient while 
analytical treatment fails and surely is more cheaper than performing experiments. 

In this work with the SOS model ||, |^, ^, ^ we would like to check how substrate temperatures 
influence qualitatively the surface morphology before the percolation limit, when only a small fraction of 
the first monolayer (ML) is deposited. 



2 Model 

The present model |[ 0, ||, |^ on a square lattice is based on random deposition (RD) enriched by a 
relaxation process of incoming at the fixed rate particles. After a random choice of the place of the particle's 
initial contact to the surface, particles migrate to one of the four nearest neighbor (NN) sites or stick on the 
place of deposition. The migration procedure is repeated L^if times, thus L^if rnay be considered a range 
of diffusion Q. In the local relaxation process, the particle tends to maximize the number of particle- 
particle lateral bonds (PPLB) similarly to the Wolf- Villain model The probability of choosing one of 
the five accessible sites is given by the Boltzmann factor exp(— iJ/fc^r), where the particle's total energy 
E = UxJx + n-yJy in all five positions depends on the number Ux (uy) of the virtually created bonds in x 
(y) direction. Jx and Jy reflect strengths of the PPLB. ks is the Boltzmann constant and T denotes the 
substrate absolute temperature. Additionally, the system tendency to increase the number of PPLB is 
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slowed down by the barrier for diffusion V: the probabihty of movement (choosing NN sites) is reduced by 
a factor exp^Vx/ksT) or exjp{Vy/kBT), depending on the considered direction of diffusion. The diffusion 
barrier V must be positive, while the negative E is compatible with the assumed tendency of the system 
to maximize the number of PPLB. 

The model control parameters were taken after Ref. (where kinetic MC simulations were per- 

formed) to mimic the realistic ratio of Vx/Vy and Jx/ Jy We assume that activation energy for diffusion 
and activation energy for breaking PPLB (in terminology of the cited paper) may correspond to the 
model control parameters V and —E, respectively. a;-direction is parallel to dimer rows, while y shows 
the direction of cross-channel diffusion on a bcc(llO) plane: 

Ag: Vx/Vy = 0.736 Jx/Jy = 9.000 
Cu: Vx/Vy = 0.793 Jx/Jy = 6.857 

3 Results of simulations 

The simulations were carried out on a 256 x 256 lattice with periodic boundary conditions. Diffusion was 
limited to Ldif = 10'^ steps. We measure the surface anisotropy in terms of the dimensionless parameter 
e= [0(1,0) -G(0,1)] /G'(0,0), where G(s) = (/i(r + s)-/i(r)) - (/i(r))2 is height-height correlation function 
[Q. The average (...) runs over all substrate sites r and the dispersion G{Q, 0) cr^ = {[h{v) — {h{v))]^) 
is defined as a square of the surface width cr, which usually describes the surface roughness. 

3.1 Influence of the substrate temperature on the surface anisotropy 

Fig. ^ shows changes of surface morphology with increasing temperature for submonolayer coverage 
9 = 0.1 ML. For the limiting case T ^ {exp{V/kBT) oo), the large diffusion barriers stop any 
migration. Thus the probability to stay in the place of deposition Pq goes to unity, while the probability 
of the migration to the NN sites Pnn goes to zero, and a film surface grown on a very co/d substrate contains 
randomly, isotropically distributed single atoms. The anisotropy parameter e tends to zero. Increase of the 
temperature changes the ratio between Boltzmann factors drastically. For low temperatures isotropically 
distributed small aggregates on the substrate appear (see the upper part of Fig. |^). In such a case 
although particles are allowed to make idif — 10^ steps, they stick after making only a few of them, 
and the surface morphology still corresponds roughly to that generated with RD rules. Then in medium 
temperatures range — see the medium part of Fig. ^ — migration in y-direction is still stopped by 
diffusion barriers, while particles are able to diffuse in a;-direction and long 1-D stripes on the substrate 
grow. The anisotropy parameter e becomes maximal. Finally, in the bottom of Fig. ^ nucleation of 
larger islands is observed with tendency of their orientation along a:-direction when the temperature is 
relatively high. For infinitely large temperature T — > cx) {exp(V/kBT) 1), the barrier for diffusion does 
not influence the particle migration. The probabilities of choosing any of the accessible sites become equal 
(Pq — Pnn — 1/5)- Particles migration is stopped only by limited to Ldif steps diffusion. Thus, for very 
hot substrates the surface morphology is similar to that generated with RD and again e = 0. The changes 
of the surface morphology may also be observed by the temperature dependence of the height-height 
correlation function G(s). In Fig. ||, as discussed earlier, the influence of the substrate temperature T on 
the surface morphology anisotropy parameter e is presented. 

3.2 Influence of the substrate temperature on the surface roughness 

For submonolayer coverage, the film contains 9L^ occupied substrate sites (with h — 1) and (1 — 6)L^ 
empty ones {h = 0). Neither infinitely large diffusion Ldif oo nor varying substrate temperature T 
can change the distribution of film heights and/or influence the single-site characteristic a — the surface 
roughness (they influence, however, two-sites characteristics such as £ — the surface anisotropy). Thus 
roughness a grows like Vd as it is predicted by Poisson process for RD case. For example, for surfaces which 
exhibit self-afRnity and, as a results of that, their roughness dynamics may be described by Family- Vicsek 
law [Q , before completing first ML — or more precisely, before substrate coverage reach the percolation 
threshold and the notion of film thickness can be applied — surface roughness dynamics corresponds to 
simple RD. Such a situation was already observed for many SOS models, and explained as an artifact of 
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Figure 1: Snapshots from simulations for coverage ^ = 0.1 ML for a 40 x 40 cut out net of the substrate. 
Dots correspond to occupied sites. Substrate temperature grows from top to bottom of figure. Columns 
correspond to Ag-like and to Cu-like parameters, respectively. 
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Figure 2: Anisotropy parameter e for different substrate temperatures T (normalized to the temperature 
To for which e is maximal) for homoepitaxial growth on (110) plane for Cu-like (x) and Ag-likc (+) model 
control parameters. Net size is 256 x 256. The substrate is covered by 6' = 0.1 ML. 



an initially flat substrate ||, |T^, |l^, |T^. For instance, in the Family model ||T^, particles after RD perform 
a local search for the column with minimal height. Thus, starting with perfectly flat substrate, and for 
submonolayer coverage, all sites are energetically equivalent and diffusion is unimportant. Note, that for 
submonolayer MBE also other law to describe dynamic scaling was proposed Q. 



4 Conclusions 

In summary, we found from computer simulations that with increasing substrate temperature a transition 
in surface morphology takes place. In contrast, varying the substrate temperature does not influence the 
film roughness for submonolayer coverage. The changes in surface morphology are due to changes of the 
relative ratio between Boltzmann factors for diffusion barriers and strengths of attraction (bonds energies) 
in two directions perpendicular to each other. 

Qualitative agreement with kinetic MC results of simulation with fully reversible aggregation SOS 
model for Ag(llO) and Cu(llO) lO, ^ |l^ and experimental observation for Cu submonolayers |l^, |2^ 
was achieved. Also STM images [|l] and hot STM movies of homoepitaxial submonolayer growth on 
Si(OOl) substrate shows similar morphology transition. For low deposition rates, 1-D islands grow well 
within a 50 K window around 530 K. At lower temperatures, mobility is decreased, and atoms deposited 
on the surface nucleate rather new islands than attaching to existing islands. At higher temperatures 
islands merge and become 2-D but remain still anisotropic . At low coverage and at room temperature 
III group elements (Ga, In) and IV group elements (Si, Ge, Sn, Pb) all form 1-D chains of ad-dimers 
|p3| . Shen and co-workers — as it was earlier stated also in the work of Mo et al. and confirmed by 
STM studies |Q — conclude that, for such a systems, not only the anisotropy in the diffusion but also 
the anisotropy in the sticking must be involved in modeling film growth. In presented here model first 
assumption is realized by different barriers for diffusion {Vx ^ Vy) while anisotropy in bonds [Jx ^ Jy) 
mimics anisotropic sticking probability. 
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